Electronic structure of self-assembled quantum dots: comparison between density 
functional theory and diffusion quantum Monte Carlo 
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We have calculated the exchange, correlation, and total electronic energy of a realistic InAs self- 
assembled quantum dot embedded in a GaAs matrix as a function of the number of electrons in 
the dot. The many-body interactions have been treated using the local spin density approximation 
(LSDA) to density functional theory (DFT) and diffusion quantum Monte Carlo (DMC), so that 
we may quantify the error introduced by LSDA. The comparison shows that the LSDA errors are 
about 1-2 meV per electron for the system considered. These errors are small enough to justify 
the use of LSDA calculations to test models of self-assembled dots against current experimental 
measurements. 
PACS Numbers: 85.30.Vw 
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I. INTRODUCTION 

The electronic structure of quantum dots resembles 
that of atoms. Both systems display three dimensional 
electronic confinement leading to level degeneracies, shell 
structure, and spin correlation to only menticw the most 
studied atomic properties of quantum dots£r0 Direct as 
well as exchange Coulomb interactions are also present in 
quantum dot systems, the importance of each depending 
on the dot dimensions .13 That is why some of the theo- 
retical methods used to calculate the electronic structure 
of quantum dots are the same as the ones used to study 
atoms and molecules. One of these methods, local spin 
density approximation (LSDA) density functional theory 
(DFT) within the effective mass approximation (EM A), 
has been widely used_to calculate the electron-electron 
interaction in dotsLTHl because of its simple implemen- 
tation and low computer demand. LSDA is an approxi- 
mate theory, well known to predict incorrectly the phys- 
ical properties of some molecules and. solids, while per- 
forming well on many other systems. lj However, experi- 
ence with LSDA on atoms and molecules may not carry 
over to quantum dots since the confining potentials and 
electron densities can be very different. As a general 
rule, density functional theory tends to work better for 
high effective density, and fails for low density systems, 
where correlation effects become important. There have 
been several investigations of LSDA and exact treatpiept 
of electron interactions in parabolic quantum dots,£2r££l 
with particular emphasis on low density systems and ef- 
fects of external magnetic fields. The electron interac- 
tions in parabolic dots have Ip^aa, treated exactly using 
diffusion Monte Car^lCMC^'Lil as well as exact diag- 
onalization methodsli3ii3. The parabolic potentials are 
popular models because they are computationally con- 
venient and, with two adjustable parameters, have had 
relative success in explaining the influence iff many-body 
effects on the electronic properties of dotsfl 



In this paper we compare EMA results obtained 
with LSDA in self-assembled InAs/GaAs quantum dots 
against an exact treatment of the EMA many-body inte*-, 
actions obtained with diffusion quantum Monte Carlo.tLZl 
These small dots, which can hold about six electrons, 
have complicated confining potentials that differ consid- 
erably from the parabolic dots considered in previous 
DMC and exact diagonalization calculations. The sin- 
gle particle states of realistically modeled pyramidal dots 
are significantly modified from the states in a parabolic 
dot, LI so many-body interactions should also differ from 
those in parabolic models. The purpose of this compar- 
ison is to quantify errors due to LSDA for calculations 
on realistically modeled InAs/GaAs dots. Although er- 
ror of LSDA is expected to be small in these systems, it 
is large enough to be a limitation to the comparison of 
model results to experiments. For the system reported 
here, we consider LSDA errors of up to 10 meV to be 
acceptable given the current precision of the models and 
experiments. Determination of whether LSDA meets this 
acceptability criterion requires careful comparison with 
many-body calculations on realistically modeled poten- 
tials. 

In these calculations we consider the ground state en- 
ergy as a function of the number of electrons in the dot, 
which is the quantity for which DFT should apply. Elec- 
tron addition and removal energies are experimentally 
measurable. These may be rigorously defined as differ- 
ences in the total energies of the ground states of the 
dot which differ by one electron. A common approach 
to calculating electron addition and removal energies is 
to use the eigenvalues of the DFT equations, which are 
well known not to be interpretable as electron addition 
and removal energies. We have addressed this issue in a 
previous papeiQ where we have shown that addition and 
removal energies are to a very gped approximation given 
by the "Slater transition rule" ,EJ which uses the eigen- 
values at one-half occupation. Thus the present work, 
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which compares total energies, is a rigorous test of the The total energy can then written as 
LSDA, and combined with our previous workQ is a test 

of the accuracy of the addition and removal energies ob- E = ^ kin + ^ cb M + Eii M + Ex M + E ° \p\ > 

tained using LSDA eigenvalues and the Slater transition w here 

rules. 

The ground state energy is the lowest eigenvalue of the ft 2 JX^ i * -i 

many-body Schrodinger equation, which we take to be ^ kin = ~~2 ^ / ( r )^(M V)ipi(r,<r)dx 



(7,2 — 1 ' 



—Wi [M- x Vi] + V cb (R) + V CC {R) \ 9{R) = EV(R)} S thc kmctic eneT SY of the Kohn-Sham orbitals, 

i=l J f 

E cb[p] = j p{r)v ch (r)dr, 



where ^(R) — ^(ri, r2, • • • rjv) is the TV-electron wave 
function, M is the electron effective mass tensor and the 
potential energy terms V c h(R) and V ee (R) are given by 

N 

V ch (R) =J2 v cb{ri) (2) 

i=i 

where w c b(r) is the potential for a single electron due to 
the offset and strain potential of the conduction band,El 
and V ee (R) is the Coulomb interaction of the conduction 
electrons with charge — e in a dielectric characterized by 
the dielectric constant e. The offset of the conduction 
band edge causes a step potential at the edge of the 
dot, which is reduced by strain, resulting in a potential 
w c b that is zero outside the dot and a roughly constant 
depth of several tenths of an eV in the interior of the dot, 
which typically has dimensions « 10 nm. Strain causes 
the electron effective mass to become anisotropic leading 
to a mass tensor given by .diag(Af ) = (to xx m yy m zz ) 
and zero off-diagonal terms .0 In the usual case of sample 
growth along the crystal direction (0 1), the electron 
masses along the plane perpendicular to the growth di- 
rection are equal, i.e. m xx = m yy . 

II. DENSITY FUNCTION THEORY AND THE 
LOCAL SPIN DENSITY APPROXIMATION 

The DFT approach to obtaining the ground state en- 
ergy is to replace the rather complex N-electron ground 
state wave function and the associated Schrodinger equa- 
tion by the much simpler ground state electron density 
p(r) and the corresponding functional forms T[p] and 
V[p] of th&,k-iaetic and potential energy operators, re- 
spectively.t^rO However, those functional forms are un- 
known, and r-approximations are necessary. In Kohn- 
Sham theoryu the density is given by a sum of densities 
of single particle orbitals, 

N 

p(r)=£E^MI 2 - (4) 

a i—1 



(•5) 



(6) 



(7) 

is the potential energy from the conduction band offset 
and strain, 

Eb\p\ =-ef p(r)0(r)dr, (8) 

is the Hartree energy, and E x [p] and E c [p] are functionals 
of the electron density. The self-consistent electrostatic 
potential </>(r) is obtained as a solution to the Poisson 
equation 

V 2 0(r) = %{v) (9) 

with the boundary condition <j>(r) — > as |r| — > oo. Mini- 
mization of the total energy defined by Eqs. (|4| - ^) yields 
a one particle Schrodinger equation which may be solved 
self-consistently with the expression for the density and 
the Poisson equation to determine the exact many-body 
ground state energy. However, the nonlocal functionals 
E x [p] and E c [p] are in general not known, so an approxi- 
mation must be made in order to proceed with the Kohn- 
Sham approach. 

In the local spin density approximation (LSDA) the ex- 
change and correlation functionals are taken to be func- 
tions of the local charge and spin density, which are then 
uniquely defined by the requirement that the approxima- 
tion be exact for the homogeneous electron gas. These 
functions have been fit to the exchange and correlation 
of the uniform electron gas.ca The validity of LSDA has 
been a subject of much research in atomic, crys tallin e, 
molecular, and parabolic quantum dot systems ,E9E!l>t3 
but has not been checked for self-assembled dot poten- 
tials, which are qualitatively like finite square wells with 
rather high electron density. The density gradients and 
shell structures of self-assembled dots cause the exchange 
and correlation energy to differ from that of a uniform 
gas. The purpose of this paper is to quantify this differ- 
ence as a function of the number of electrons in a realis- 
tically modeled dot. 

III. QUANTUM MONTE CARLO METHODS 

The simplest QMC method is variational Monte Carlo 
(VMC), in which expectation values for a trial wave func- 
tion are evaluated using a Metropolis algorithm. Using 
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the variational principle of the ground state energy, pa- 
rameters in the wave function can be optimized to min- 
imize the total energy, thus providing an approximation 
to the ground state energy and wave function. In our cal- 
culation we use VMC to optimize the guiding wave func- 
tion for our diffusion Monte Carlo algorithm, described 
below. 

We also use Monte Carlo integration to evaluate the 
exact exchange energy. In the Kohn-Sham formalism, 
the exchange energy is defined as the difference 



E x [n] = (*ks|Voc|*ks) -E u [p], 



(10) 



where |^ks) is a Slater determinant of the true Kohn- 
Sham orbitals; those that give the exact density. We 
assume that LSDA is accurate enough to provide approx- 
imations to the Kohn-Sham orbitals, i. e. that the LSDA 
density is close to the exact density. We then evaluate 
the integral (^ksIKcI^ks) using a Slater determinant of 
our calculated LSDA wave functions. 

Diffusion Monte Carlo is a stochastic method that is 
able to project the N-electron ground state wave func- 
tion <I>o from a trial wave function Wt- A position basis 
is used to described the state of the system. A good trial 
wave function is important because the variance of the 
Monte Carlo sampling decreases as the trial wave func- 
tion approaches the true ground state wave function. We 
chose the trial wave function to be a product of Slater 
determinants for spin up and spin down electrons with a 
Jastrow factor that introduces correlation, 



* T (R) = det[0 k (rO] x exp[- ^u(^)]. 



(11) 



KJ 



For convenience, we take the single particle wave func- 
tions #k(r) to be the non- interacting eigenstates for the 
external potential, and choose the Jastrow factor u(r) — 
—ar/(l + br) with a — m/e and b = 1/2. 

Following the method described in Ref. [ |Ti| , we 
used importance sampling from a function $(i?) = 
*S?t(R)&o(R)- The anisotropic mass is taken into account 
using the mass tensor, M. The projection of a state R to 
the ground state distribution ^t^o is accomplished by 
repeatedly applying the projection operator e~ Hr , each 
time sampling a new configuration. Here H is the Hamil- 
tonian and t is a parameter chosen to be small so that 
the projection operator is approximated bji!3 



(R\^ T {R)e- HT ^\R')\R') 
= exp{-[E L (R) - E ]t} 
[(2 7 rr) 3 detM- 1 ]t 



(12) 



J|cxp 



(jCi -n- TV d;i ) ■ M ■ (i-j - n - rv d - 
2r 



where Eq is a constant parameter to control the popula- 
tion, E L (R) = •fy T (R)~ 1 H 1 S> T (R) is the local energy, and 
v d = M _1 • Vlog^T- The algorithm is thus: (1) Start 
with an ensemble of configurations (walkers) distributed 



GaAs 




FIG. 1. Schematic representation of the self-assembled dot 
structure used in the present work. 

by |^t| 2 , (2) propagate each configuration with a drift 
Vd and Gaussian displacement with covariance matrix 
(Ty = t, (3) reweight each configuration by a fac- 
tor exp{ — (£l — E )t}, (4) repeat, collecting statistics 
once the steady state distribution is reached. The exact 
ground state energy is obtained by sampling the local 
energy. We use branchings at each step to improve the 
efficiency of the process. 

A complication known as the fermion sign problem 
arises when applying diffusion Monte Carlo to fermions. 
Electron exchange introduces negative signs into the pro- 
jection operator, which decreases the efficiency of the 
Monte Carlo sampling. In the present discussion the 
short time approximation, Eq. (|12|), breaks down when 
walkers cross nodes of 'J't- We handle this problem with 
the fixed node approximation,!^ in which walkers are 
given a weight of zero when they cross a node of \&t- 
This approximation has a lower bound property, so the 
true ground state energy can only be lower than the fixed 
node energy, and the exact choice of \&t = would give 
the exact Fermionic ground state energy. 

The correlation energy can be deduced from a knowl- 
edge of all other energies, 

E c [p] = E iot — Skin — E c b — En — E x , (13) 

or equivalently 

E c [p] = (* |Vco|*0> - (*Ks|Vee|*KS>. (14) 



IV. NUMERICAL COMPARISON 

Using DFT and DMC we have calculated the total 
ground state energy of the system of electrons, defined 
for DFT by Eq. |, and defined within DMC as the fixed 
node ground state energy, as defined in the previous sec- 
tion. Figure [l] shows the system we have considered, a 
pyramidal shaped InAs quantum dot in a GaAs matrix. 
However, our method of analysis is not restricted to this 
particular shape, and we expect our results to be valid 
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for truncated pyramidal or lens shaped dots. We have in- 
troduced several simplifications in the model. Since our 
DMC scheme does not allow for a position dependent di- 
electric constant, we have used only one value, e = 14eo, 
for both GaAs and for InAs. In the absence of strain the 
GaAs - InAs conduction band offset between the two ma- 
terials is 770 meV. To partially take into account strain 
effects, we have assumed a constant shift of 400 meV in 
the InAs conduction band, and no strain effect in the 
GaAs conduction band, resulting in a conduction band 
offset of 770meV - 400meV = 370meV. We also take the 
electron effective masses to be constant in each material. 

In order to study a realistic potential, we have used a 
nonuniform grid basis for the LSDA calculation and the 
solution of the Poisson equation. For the Schrodinger 
equation we set the wave functions to zero at the edge 
of the grid, which is reasonable due to the exponential 
decay of bound states, but the Poisson equation requires 
more care. We used a multipole expansion up to the 
quadrupole term to simulate the boundary conditions at 
infinity. This difficulty with the boundary conditions of 
the Poisson equation is one source of error in our com- 
parison. We used cubic interpolation to map the gridded 
wave functions and conduction band potential to continu- 
ous functions for Monte Carlo evaluation of the exchange 
energy, for determination of the trial nodes in the DMC 
calculation, and for guiding wave functions in the DMC 
calculation. Although the cubic interpolation of the po- 
tential allows us to compare LSDA and DMC for similar 
external potentials w c b, errors in the kinetic energy op- 
erator are more difficult. The use of a finite difference 
approximation to the Laplacian in the computation of 
the LSDA solution creates an operator that has no sim- 
ple counterpart in a continuous formulation of the prob- 
lem. The error between the physical Laplacian operator 
and the artificial finite difference operator can be made 
arbitrarily small by the use of finer grids, but computer 
memory and CPU time constraints caused this to be a 
significant source of error in our comparison. The size of 
the errors are several meV, which are comparable to the 
rather small LSDA errors, especially in the correlation 
energy. Below we discuss a way we estimated this error 
so that we can compare correlation energy. 

For the set of calculations we describe now, we have 
taken all masses to have the isotropic value of 0.05m e , for 
both the InAs and GaAs regions. The case of anisotropic 
masses will be discussed later. We have also accounted 
for numerical discrepancies between the two calculations 
due to grid interpolation error. To estimate this error we 
have compared eigenvalues of single electron states with 
both methods, and found that the DMC single particle 
eigenvalues typically lie about 0.25 meV above the eigen- 
values computed by the grid method used in the DFT 
calculation. We attribute that difference to a systematic 
error in the finite difference kinetic energy operator and 
shift up our LSDA calculations to compensate for this. 

Figure ^| shows the values of the different components 
of the total energy as a function of dot occupation. Dif- 
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FIG. 2. Contributions to the energy of a 200 x 100 A pyra- 
midal dot as a function of occupation TV, calculated by both 
LSDA and QMC. Differences between LSDA and QMC are 
not apparent on the energy scale of this figure. This clearly 
shows the small effect of interactions beyond the Hartree 
(mean field) level. 



ferences between LSDA and QMC are not apparent on 
the scale of the figure. This figure clearly shows that the 
external potential energy and kinetic energy are much 
larger than the interaction energies. In other words the 
effects of interactions enter as a perturbative correction 
to the non-interacting system. The reason for this can 
be seen from the scaled electron density. The energy and 
length scales for the electron interaction are scaled by the 
dielectric constant and mass, so that the effective Bohr 
radius is a*, = e/m*ao w 150A and the effective Hartree 
is Ha* = e~ 2 m*Ha s=s 7meV. If we estimate the density 
by assuming the N electrons uniformly occupy the inte- 
rior of the dot, we find an effective conduction electron 
density of r s w OAdN" 1 ^ 3 . To see the consequences of 
such a high effective density, consider the uniform elec- 
tron gas at r s = 0.6. The electxon gas has a ground state 
energy expansion for small r s ,E3 

E = 2.20997-7 2 - O.gieS^T 1 - 0.094 + 0.0622 ln(r s ) + • • • , 

(15) 

where the first term is the kinetic energy, the next term 
is the exchange, and remaining terms are correlation en- 
ergy. For the case of six electrons in the dot r s w 0.25 and 
the expansion gives -Ekin ~ 1440meV, E x w 150meV, and 
E c w 0.6meV. Although the comparison between these 
very different electronic systems cannot be pushed too 
far, this does show that our exchange and correlation en- 
ergies are reasonable for this electron density. The lead- 
ing effect of the interactions is the Hartree energy, with 
small corrections for exchange and very small correlation 
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FIG. 3. Exchange energy of a 200 x 100A pyramidal dot 
as a function of occupation N, as determined by LSDA and 
QMC. Also shown is the self-energy correction, which makes 
a large contribution to the exchange energy. 



corrections. 

The most important contribution of LSDA in this sys- 
tem is the exchange term, and we plot the comparison of 
the LSDA exchange energy in Figure [|. For finite sys- 
tems, the definition for exchange includes.the correction 
for self-interaction in the Hartree energyj^j which we also 
show in the figure. Merely correcting for self-interaction 
will recover at least 75% of the exchange energy, and 
LSDA is able to do better, recovering about 90% for an 
error in the exchange energy of less than 5 meV. As is 
well known in many other systemsE3, the error in the lo- 
cal approximation to the exchange is the largest error in 
LSDA. 

By comparison the correlation energy is much smaller 
and its accuracy is difficult to assess due to grid errors. 
Correlation is defined as the difference in total energies, 
Eq. ([l4|), and in this case the total energy difference is 
less than 1% . Thus our estimates of the exact corre- 
lation energy are fairly uncertain. In Figure |] we show 
our best attempt at checking the LSDA correlation in 
the quantum dot. Since the LSDA wave functions are 
not the true Kohn-Sham orbitals (due mostly to grid er- 
rors) the correlation is overestimated due to relaxation 
of our interpolated LSDA states. As stated earlier, we 
correct for this by performing single particle calculations 
to estimate the relaxation of our LSDA states, and have 
removed this contribution from our calculated correlation 
energy. These results are shown as "QMC" and "QMC 
with corrections" in the figure. We find that LSDA gives 
a reasonable estimate of the magnitude of the correlation 
energy in this system. 

Our comparison of the energy errors from the LSDA 
are summarized in Figure |. The largest error is due to 
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FIG. 4. Correlation energy of a 200 x 100A pyramidal dot 
as a function of occupation N, as determined by LSDA, QMC, 
and QMC with corrections for grid interpolation error, as ex- 
plained in the text. LSDA gives a reasonable approximation 
to this small quantity. 

LSDA underestimating the magnitude of the exchange 
energy, E x , which is a 6 meV error for six electrons in the 
dot. There is a small error in the correlation energy, E Cl 
which slightly compensates for the error in the exchange. 
The error in the Hartree energy i?Ha is due primarily 
to truncation of the multipolc expansion at quadrupole 
terms for determination of the boundary conditions of the 
Poison equation. There are also errors in -Ekin and E c ^ 
(not shown in the figure) due to the grid, in particular 
the discretization of the kinetic energy operator discussed 
earlier. We therefore find that LSDA leads to an error in 
the total energy of up to 6 meV, which in our calculation 
is partially compensated for by grid errors in other energy 
terms, leading to an overall error in the total energy of 
no more than 2 meV. 

Because of strain, the electron effective mass in the 
plane parallel to the base of the djt is usually different 
from the mass perpendicular to it .El The LSDA does not 
explicitly account for mass anisotropy in the exchange 
and correlation functionals, thus for the calculation of 
the exchange-correlation potential V xc we have assumed 
a constant average mass given by 

m xc = (mjjmjj 1/3 , 



(16) 

where mil is the in-plane mass, mn = m xx = m yy , and 
m± is the perpendicular mass, m± — m zz . The choice 
of the average mass m xc is somewhat arbitrary. We have 
performed a test of this approximation by calculating the 
exchange and correlation energy for an anisotropic sys- 
tem using (i) LSDA with m xc as given above, (ii) LSDA 
with m xc optically averaged, m^. 1 = (2/m\\ + l/m±)/'S, 
and (iii) DMC, which can explicitly treat anisotropic 
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FIG. 5. Difference in energies between LSDA and QMC 
(-Elsda — -Bdmc) for a 200 x 100A pyramidal dot as a function 
of occupation N: (a) error in the exchange energy, £ x , due 
mostly to LSDA, (b) error in correlation energy, E c , due to 
grid interpolation error and LSDA, with some grid corrections 
as described in the text, (c) error in Hartree energy Ehb,, due 
to grid interpolation error and the multipole expansion, and 
(d) error in the total energy. The values in (a), (b) and (c) do 
not add up to the value in (d) due to errors in E c b and Skin 
(not shown). 

mass. Figure ^ shows the calculated exchange and corre- 
lation energies for several electron mass ratios, which we 
have again fixed to constant across the InAs and GaAs. 
We find that either choice is acceptable, as both give er- 
rors comparable to the isotropic mass case. 

V. CONCLUSIONS 

In conclusion, we have shown that DFT offers an accu- 
rate approximation for the ground state energy including 
many-body interactions in small self-assembled quantum 
dots, while providing a simple and fast means of model- 
ing systems containing several electrons. This is in large 
part due to the small size of the self-assembled quan- 
tum dots in relation to the effective Bohr radius. In this 
regime the largest error is in the local approximation to 
the exchange, and we have verified that the errors due to 
this approximation are small. This is quantified in Fig- 
ures |^ and [| where we see that the exchange energy for 
six electrons in the dot is 72 meV with an LSDA error 
of 6 meV, compared to a total energy of 1490 meV. We 
expect this applicability of LSDA to hold in general for 
small quantum dots of various shapes, but we emphasize 
that LSDA is expected to become progressively worse if 
the dot becomes much larger than the effective Bohr ra- 
dius. Also, the use of the EMA Schrodinger equation, 




1 1.5 2 

m x/ m z 

FIG. 6. Values of the (a) exchange and (b) correlation 
energy as a function of mass anisotropy for six electrons in 
a 200 x A pyramid dot. Triangles are DMC results, x's are 
LSDA results with — (2mjj _1 +m]| 1 )/3, and squares are 

LSDA results with m xc = (m| * mi) 1 ' 3 . 

Eq. (Q), although reasonable for the .chosen problem, 
is predicted to fail for very small dots.ES Although our 
model assumes pure InAs dots surrounded by GaAs, in. 
some systems indium and gallium may in fact intermix.EZI 
Again, we expect the relative magnetude of exchange and 
correlation energies and the acceptability of LSDA to be 
unchanged by intermixing or alloying of the dots. 

Based on this comparison we conclude that the ad- 
dition and removal energies found by EMA-LSDA cal- 
culations in realistic dotsQ are accurate to 1-2 meV per 
electron as far as the many-body interactions within the 
EMA are concerned. Thus, comparison of such calcu- 
lated charging energies with experiment can be consid- 
ered a direct test of the models and the uncertainties 
in the analysis of experiments.^ Furthermore, we have 
carefully tested the grid errors and concluded that ener- 
gies are accurate to « ±10 meV including all errors. 

We would like to point out that application of LSDA to 
a system of several coupled self-assembled dots may have 
difficulties. Systems that have weak intersite transitions 
along with intrasite interactions are a well-known many- 
body problem where LSDA is known to fail. The prob- 
lem is closely related to the H2 molecule, where LSDA 
describes well the electronic energy at equilibrium dis- 
tance, but gives an incorrect broken symmetry solution 
if the atoms are pulled apart beyond a critical distance. 
The problem of weakly coupled dots would be of partic- 
ular importance for future studies. 
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